Screened and Unscreened Phases in Sedimenting Suspensions 
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A coarse-grained stochastic hydrodynamical description of velocity and concentration fluctuations 
in steadily sedimenting suspensions is constructed, and analyzed using self-consistent and renormal- 
ization group methods. We find that there exists a dynamical, non-equilibrium phase transition 
from a "unscreened" phase in which we recover the Calflisch-Luke (R.E. Calfiisch and J.H.C. Luke, 
Phys. Fluids 28, 759 (1985)) divergence of the velocity variance to a "screened" phase where the 
velocity fluctuations have a finite correlation length growing as <^ _1 ' /3 where <j> is the particle volume 
fraction, in agreement with Segre et al. (Phys. Rev. Lett. 79, 2574 (1997)) and the velocity variance 
is independent of system size. Detailed predictions are made for the correlation function in both 
phases and at the transition. 
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Sedimentation Q is a rich and complex phenomenon 
in suspension science and a frontier problem in nonequi- 
librium statistical mechanics. The average sedimenta- 
tion speed i> S cd of solute particles drifting down in a sol- 
vent is determined by balancing the driving force (grav- 
ity) against the dissipative force (viscous drag). Giant 
non-thermal fluctuations in the velocity and concentra- 
tion fields in a steadily settling suspension, observed even 
for non-Brownian systems, have been a puzzle for some 
years. Cafiisch and Luke (CL) ||] showed, for steady sed- 
imentation in a container of smallest linear dimension L, 
that the assumption of purely random local concentration 
fluctuations led to velocity fluctuations with a variance 
(v 2 ) ~ L. Most experiments, however, find no depen- 
dence of (v 2 ) on L although Ladd's simulations ||] 
and the data of Tory et al. Q appear to be consistent 
with CL. 

In this Letter we propose a resolution of this puzzle by 
means of a set of coarse-grained, fluctuating nonlinear hy- 
drodynamic equations for the long-wavelength dynamics 
of concentration and velocity fluctuations in a suspension 
settling steadily in the — z direction, at vanishingly small 
Reynolds number. Our theory is similar in spirit to the 
Koch-Shaqfeh (KS) || "Debye-like" screening approach 
but differs in several important details and predictions. 

The central conclusion of our study is that there are 
two qualitatively distinct nonequilibrium phases for a 
sedimenting suspension. In the "unscreened" phase (v 2 ) 
diverges as L, as in CL and, in addition, concentration 
fluctuations with wavevector k = (kj_, k z ) relax at a rate 
cx fc 1 / 2 . The "screened" phase is characterized by a cor- 
relation length £ similar to that predicted by KS such 
that (t; 2 ) ~ L for L <C £ and (v 2 ) ~ £ for L > £. Deep 
in the screened phase we predict £ ~ where <f> is 

the particle volume fraction. This is in agreement with 
the experiments of Segre et al. || , but not with KS || . 
The relaxation rate in the screened phase is independent 



of k for k z — and k^ — > 0. Detailed, experimentally 
testable expressions for the structure factor and velocity 
correlations in the screened phase are presented after we 
outline our calculations. The two phases are separated 
in our "phase-diagram" (Fig. 1) by a striking continuous 
nonequilibrium phase transition where £ diverges at least 
as rapidly as (K — K c ) 1 as a control parameter K is 
decreased towards a critical value K c . 

The hydrodynamic equations we used to arrive at these 
results are 

Br 

— + v • Vc = [£> jVi + D 0z V 2 ]c + V • f (r, t) (1) 
and 

?7V 2 ^(r,i) = m R gP iz c{r,t), (2) 

where c(r, t) and v(r, t) are the fluctuations about the 
mean concentration cq and the mean sedimentation ve- 
locity —v sc dZ respectively. We justify these equations 
briefly below; for a more detailed discussion we refer 
the reader to Ref. ||. Eq. |l] is the anisotropic ran- 
domly forced advection-diffusion equation with bare uni- 
axial diffusivities (Dq z , Dq±) and a random stirring force 
f(r,i) 0]. The Stokes equation, Eq. |, which expresses 
the balance between the driving by gravity and the dis- 
sipation by the viscosity 77, describes how the concen- 
tration fluctuations produce velocity fluctuations. Here 
mug is the buoyancy-reduced weight of a particle, while 
the pressure field has been eliminated by imposing in- 
compressibility via the transverse projection operator 

Hydrodynamic equations such as Eqs. [I] and || arise 
from a coarse-graining of the microscopic equations of 
motion. The latter, for the main case of interest here, 
viz., non-Brownian suspensions at zero Reynolds num- 
ber, are the deterministic equations of Stokesian dynam- 
ics for N hydrodynamically coupled particles, and are 
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known to be chaotic pij . The noise, or random stirring 
current f(r, t) and the diffusivities in Eq. [I] represent a 
phenomenological description of the deterministic chaos 
at length scales below the coarse-graining length £ (which 
must be large compared to the particle radius a). We 
use these hydrodynamic equations to predict the veloc- 
ity and concentration fluctuations at length scales large 
compared to £ driven by the random stirring at short 
distances. 

We assume, as is reasonable, that f(r, t) is Gaussian 
white noise with uniaxial symmetry: 



(/i(r J f)/ J -(r' l f / )) = 2c N* j S(r - r')S(t - t') 



(3) 



with an anisotropic noise amplitude Nq = Nq±S^ + 
N 0z S^, where and S^j are the projectors along 
and normal to the z axis, respectively. Because of 
the nonequilibrium origin of the noise and diffusion 
constants, we may not jl2| assume that Nq±/No z = 
Dqj_/Dq z as would be true for the Langevin equation 
of a dilute suspension at thermal equilibrium. Note that 
no correlations have been fed in via the noise: any that 
emerge in the long-wavelength properties are a result of 
the interplay of advection and diffusion. 

Let us now consider the nature of the spatio-temporal 
correlations implied by Eqs. [l] and ^|. We will focus on 
the structure factor for concentration fluctuations 



S(q) 
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from which the velocity structure factor can be derived 
through Eq. || If we ignore the advective nonlinearity 
v ■ Vc, then S'(q) can be computed by straightforward 
Fourier transformation of Eq. |l|, resulting in 



5(q) = So(q) 



N 0± qj + N 0z q 2 z 
D 0± q 2 ± + D 0z q 2 z 



(5) 



Using Eq. || in Eq. || we can compute (v 2 ) as a function 
of the system size L with the result: 



(v 2 ) 



S(q) 



L, 



q>l/L 1 



(6) 



In other words, neglecting large-scale advection by the 
velocity fluctuations leads to the CL result. 

To include the effect of the advective nonlinearity we 
have performed a self-consistent mode coupling calcula- 
tion fTJ on Eqs. Our results can be expressed in 
terms of a renormalized relaxation rate 



R(q) = D ± (q)q 2 ± + D z (q)q 2 z + T(q) 

and a renormalized structure factor of the form 
N ± (q)q 2 ± + N z (q)q 2 



S(q) = 



R(q) 



(7) 



(8) 



The quantities D z .±(q) and N Zt ±(q) represent renormal- 
ized diffusivities and noise amplitudes Jl3| . But, most 
importantly, the advective nonlinearity to lowest-order 
perturbation theory leads to an additional term in the 
renormalization of the relaxation rate which is of the 
form r(q) = ^{q)q\/ q 2 ■ Starting from the stochastic 
hydrodynamic equations, Eqs.|l|-||, it turns out that the 
amplitude of this singular contribution becomes a con- 
stant, lim g ^o7(q) oc I(Pn,Pd), which depends on the 
anisotropy ratios of the noise and diffusivity coefficients 



a Nx , R D± 



(9) 



In particular I(/?/v, (3d) is proportional to /3/v — Pd- and 
consequently may change sign upon varying the noise and 
diffusivity ratios. For I([3n : I3d) < this would lead to 
exponentially growing concentration fluctuations in the 
limit of long wavelength. Here we do not pursue this 
intriguing possibility further but instead restrict our at- 
tention to I(PN,Prj) > 0, for which the model can either 
be treated within dynamic renormalization group theory 
or using self-consistency methods. 

We start our discussion at the borderline of stability, 
Pn = Pd- For these parameter values it can be shown 
that the fluctuating hydrodynamic equations describe a 
dynamics which obeys detailed balance j^]: the advec- 
tive nonlinearity does not affect the equal-time correla- 
tions, and S(q) in particular is just the constant N±/Dj_. 
There are singularities in N± jZ and Dx,z which we dis- 
cuss later. 

For 0n > Pd, detailed balance is violated and a sin- 
gular diffusion term T(q) is generated within perturba- 
tion theory. In order to analyze the dynamics in this 
regime we use one-loop self-consistent theory (mode cou- 
pling theory) and arrive at the expression 



r(q) = c 



m R g 
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S(k)] 



k 2 q 2 



R(k) + R(q - k) 



(10) 



with R(q) given by (0), and similar self-consistent in- 
tegral equations for D±(q), D z (q), N±(q), and N z (q). 
We find that there are two types of iteratively stable so- 
lutions to these coupled self-consistent equations: those 
with 7(5 — > 0) > 0, which we obtain below the solid line 
in the phase diagram spanned by the two anisotropy ra- 
tios ("screened" phase in Fig. 1), and those with 7(5 — + 
0) = 0, which arises for values of the anisotropy param- 
eters that lie above the solid line and below the dashed 
line of the same figure, i.e., in the "unscreened" phase. 
Note that within the self-consistent theory the line in the 
phase diagram where 7(5 = 0) changes sign (solid line) 
has shifted with respect to the result of the one-loop per- 
turbation theory discussed above (dashed line). 
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Screened Phase: In the screened phase, T(q) is of the 
form jq^/q 2 in the small q limit, with 7 a finite constant. 
This implies that the structure factor at small wavenum- 
ber becomes 



S(q) 



N 



D 



(11) 



with iVj_ jZ and D±, z constants. From Eq. [Tl]we can define 
a correlation length £ = (D^/j) 1 ^ 2 such for q±_ 3> l/£ 
the structure factor is not significantly affected by ad- 
vection. On the other hand, for q± <C l/£ the in-plane 
structure factor reads S(q±,q z = 0) ~ (N±/j) q\, while 
S(q± = 0, q z ) ~ (N z /D z ). Physically, this means that at 
long wavelength advection strongly suppresses in-plane 
concentration fluctuations. 

Using Eq. [ll] in conjunction with Eq. ^, one finds that 
for length scales L less than £, (v 2 ) cx L, consistent with 
CL, while for L large compared to £, (v 2 ) oc £. Veloc- 
ity fluctuations on length scales small compared to £ are 
thus highly correlated while they become uncorrelated at 
larger length scales. 

Deep inside the screened phase, i.e., for large 7, the 
renormalization of the diffusion and noise parameters is 
negligible and we can explicitly compute 7, and thus £, 
by inserting Eq. || in Eq. [h] using the bare values for the 
N's and D's. We find 



. m R g x 
1 T)D ' 



-2/3-1/3 

c 



2 



-1/3 



(12) 



where for simplicity we have set Dq± = D$ z = D. Ac- 
cording to Eq. [l2|, the correlation length increases as we 
decrease the /?jv parameter (which could be done by in- 
creasing the thermal noise amplitude) and diverges at 
(3n = 2. Strictly speaking, as (3n ~> 2, the diffusiv- 
ity corrections are no longer negligible, and the actual 
divergence of £ is probably stronger than (|l2]), and oc- 
curs at a larger value of (3m- An explicit analytical (but 
lengthy) result for the correlation length £ can also be 
obtained throughout the screened phase as a function of 
both anisotropy parameters || and the phase boundary 
can also be computed. The phase boundary resulting 
from this result is shown in Fig. 1 as the solid line sep- 
arating the screened from the unscreened phases. The 
dashed line in the figure corresponds to the set of parame- 
ter values where the hydrodynamic equations correspond 
to a Langevin dynamics in thermal equilibrium. 

Unscreened Phase: As already noted above, the hy- 
drodynamic equations obey detailed balance JT5| along 
the line (3n — Pd in the phase diagram. As a conse- 
quence the ratio of noise to diffusivity can be identified 
as a direction-independent "noise-temperature" . Fur- 
thermore, the structure factor S(q) becomes a constant 
Dj_/Nj_ and we recover the CL result. In conjunction 
with an exponent identity resulting from Galilean invari- 
ance this is enough to determine the dynamic exponent 




FIG. 1. Dynamical phase diagram for sedimentation. Be- 
low the solid line in the parameter space spanned by the 
anisotropy factors for the noise and diffusivities velocity and 
concentration fluctuations have a finite screening length in the 
limit vanishing wavevector. This region is called "screened" 
above. In the upper region called "unscreened" the screening 
length becomes infinite. The dashed line represents the set 
of values for the anisotropy factors where the hydrodynamics 
obeys detailed balance. The inset shows the behavior of the 
phase boundary in the limit of large noise and diffusivity in 
the vertical direction as compared to the horizontal plane. 



exactly, z — d/2 — 1. This implies that the diffusiv- 
ities and noise amplitudes scale as q~ t l 2 = q^l 2 for 
long wavelength. Even though there are now singular 
corrections to D z .j_(q) and N z .±(q), the anomalous T(q) 
term is zero. For parameter values in the regime between 
the dashed line (detailed balance line) and the solid line, 
which marks the location of the nonequilibrium phase 
transition, renormalization group methods may be used 
to determine the renormalization of the noise and diffu- 
sivity amplitudes. In view of the results from the above 
self-consistency calculation (7 = in the unscreened 
phase) and the exact results at the detailed balance line 
it is quite likely that the resulting renormalization group 
flow will tend towards a fixed point which obeys detailed 
balance. We leave the details of such an investigation for 
a future publication ||. 

The analysis of our hydrodynamic equations thus con- 
firms that screening can suppress the CL divergence of 
(v 2 ) with L, as argued by KS, while it allows for a sec- 
ond, unscreened phase. This result may help explain 
the conflicting results on (v 2 ) obtained by different work- 
ers [§-0. The self-consistent structure factor, Eq. || wc 
obtained differs significantly from the one proposed by 
KS. Experimental test will thus be of considerable im- 
portance. Measurements of S(q), for example by PIV 
]l6t (Particle Imaging Velocimetry) , would constitute the 
most direct test of the theory since our prediction that 
S(<i±, q z — 0) oc q\ does not hold in the KS description. 
Detailed measurements of S(q) for sedimenting solutions 
are not yet available. However, Segre et al. || do report 
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that the size-dependence of the amplitude (v 2 ) of the 
velocity fluctuations depends on a characteristic length 
scale £s such that (v 2 ) oc £s for length scales L ^> f$ 
while for L -C £5, (v 2 ) grows with L. They report that 
£s ~ a</) -1 / 3 with <f> the particle volume fraction. 

Our correlation length £, in Eq. [l^, has the same phys- 
ical interpretation as £g. Deep in the screened phase, i.e., 
for /(/3jv,/3d) ^ 0, c; can be written as: 

£(0) ~ (mia/vDyWat-Wlip^PD)- 1 / 3 (13) 

On scaling grounds, we expect that D oc (5u RM sC with 
S v RMS the root mean square of the velocity field fluctua- 
tions. Experimentally, S v RM s£ is found to be independent 
of volume fraction 0. In that case, Eq. [l3| reproduces the 
experimentally observed volume-fraction dependence, in 
contrast to KS g. It should be noted that this volume 
fraction dependence of the correlation length implies that 
there is a fixed number of colloids within a correlation 
volume independent of volume fraction. 

The observation of a transition from the screened to 
the unscreened phase would obviously be the most con- 
clusive evidence supporting our theory, in particular if 
the transition were accompanied by a divergence of the 
velocity fluctuation correlation length. Even in the ab- 
sence of such direct evidence, the observation of screened 
behavior combined with our theory requires that the 
anisotropics in the noise and diffusivity lie in the lower 
region of our dynamical phase diagram, Fig. 1. A com- 
plete test of our theory thus requires measurement of the 
N and D parameters. These could be obtained from the 
measurement of the steady-state static structure factor 
S'(q), e.g. by particle imaging or light scattering experi- 
ments both along the z direction and in the x — y plane, 
coupled with tracer diffusion measurements. 

Finally, it would be interesting to vary the effective 
noise and diffusion constants in a controlled manner in 
an experiment. While there is, as yet, no method to cal- 
culate these constants directly from a microscopic theory 
it is reasonable to expect that by decreasing the Peclet 
number (i.e., increasing the role of isotropic thermal dif- 
fusion) one could drive the sedimenting system into the 
unscreened phase. Thus by repeating the experiments of 
Segre et al. || with colloids that are more nearly density 
matched to the solvent one could test our prediction of a 
transition to an unscreened phase. 
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